Exploring local quantum many-body relaxation by atoms in optical superlattices 
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We establish a setting — atoms in optical superlattices with period 2 — in which one can experimentally probe 
signatures of the process of local relaxation and apparent thermalization in non-equilibrium dynamics without 
the need of addressing single sites. This opens up a way to explore the convergence of subsystems to maximum 
entropy states in quenched quantum many-body systems with present technology. Remarkably, the emergence of 
thermal states does not follow from a coupling to an environment, but is a result of the complex non-equilibrium 
dynamics in closed systems. We explore ways of measuring the relevant signatures of thermalization in this 
analogue quantum simulation of a relaxation process, exploiting the possibilities offered by optical superlattices. 



Is it possible to consider the relaxation of a closed quan- 
tum system to an apparently equilibrated state? In contrast to 
the deep understanding we have of equilibrium quantum sta- 
tistical mechanics, non-equilibrium relaxation processes are 
far from being fully understood. Specifically, the question of 
how Gibbs or relaxed states emerge dynamically is a question 
that is receiving a lot of attention recently [1-9]. Part of the 
reason for the renaissance in the study of questions of non- 
equilibrium dynamics of quantum many-body systems stems 
from the fact that systems have become available that promise 
to make such issues amenable to experiment [10]: The states 
of cold atoms in optical lattices can be manipulated with a 
high degree of control, offering a testbed for questions of non- 
equilibrium dynamics. 

The setting of interest in this work is the one of a sudden 
quench [1-9, 11, 12]. Starting in the ground state of a lo- 
cal many-body Hamiltonian, system parameters are suddenly 
changed such that the old state is no longer an eigenstate of 
the new Hamiltonian, generating a non-equilibrium situation. 
The dynamics of the system is then monitored in time. It has 
been conjectured that in such a non-equilibrium situation, one 
may — in some sense — arrive at the maximum entropy state 
consistent with the expectation values of the constants of mo- 
tion fixed by the initial state @], also referred to as a general- 
ized Gibbs ensemble fill [Till . This is appealing as it parallels 
Jaynes' approach to equilibrium statistical mechanics. Yet, of 
course, if the system can be meaningfully treated as a closed 
quantum system, one cannot expect the entire system to relax, 
as an initially pure state will remain so in time [2-4] . The en- 
tire information of the initial condition is still present in the 
system, albeit in a very dilute fashion. 

This is, however, by no means inconsistent with the expec- 
tation that the system may locally appear to be relaxed [2-5]. 
Locally, such a relaxation can well be true: Any subsystem 
may appear to be in a maximum entropy state under the con- 
straints dictated by the constants of motion, and remain so for 
an arbitrary amount of time. Ref. J3] introduces an instance 
in which this local relaxation of subsystems, referred to as 
the local relaxation conjecture, can actually be proven rigor- 
ously to hold: When quenching a state in a Mott phase in a 
Bose-Hubbard system to the non-interacting deep superfluid 
phase, the reduced state of consecutive sites converges (in 
trace-norm) to a maximum entropy state consistent with the 
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FIG. 1: Sketch of a local relaxation starting from an initial condition 
of bosonic atoms being present or absent, in even and odd sites of 
a one-dimensional Bose-Hubbard system, achieved by imposing a 
superlattice to an optical lattice. 

constants of motion J2j] — without having to invoke a time av- 
erage. The intuition is that the quench creates local excitations 
that travel through the lattice at a finite speed ESQ. Their 
incommensurate influence then leads to a relaxation without 
environment. This intuition is corroborated in Ref. 0], where 
local relaxation for Gaussian initial states is shown. 

Unfortunately, while many situations that give rise to local 
relaxation may be generated, a challenge so far unresolved is 
to actually probe signatures of local relaxation: Demonstrat- 
ing local relaxation appears to necessitate local addressing, a 
requirement that poses a great experimental challenge in sys- 
tems of atoms in optical lattices. 

This dilemma will be resolved in this work, employing a 
simple yet promising idea: We will make use of a periodic 
setting; but we are by no means obliged to stick to period 1: 
The idea is that we can well make use of a period-2 setting, 
exploiting optical superlattices. This approach will open up a 
way to quantitatively explore local relaxation effects in exper- 
iment, without the need of addressing single sites at any point. 
The period-2 will allow for observing most of the relevant 
signatures. To demonstrate the validity of this idea, we will 
make use of analytical as well as numerical methods, based 
on a time-dependent density-matrix renormalization-group (t- 
DMRG) approach. Our findings provide a simple guideline to 
what is to be expected in realistic experimental situations. 

Proposed experimental setup. - Ultracold atoms in optical 
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lattices provide a great deal of control over the system's pa- 
rameters ifioll . In particular, sudden quenches of parameters 
are accessible, and on experimental timescales, systems can 
be treated as essentially closed. To address the local detection 
problem, we prop ose to study local relaxation using optical 
superlattices II 1611 . We follow the setup recently realized by 
Bloch and coworkers [fl6tl . considering bosonic 87 Rb atoms 
in a period-2 optical superlattice geometry. This experimen- 
tal setting allows for changing the relative intensity of the two 
optical lattices, shifting their relative position, and coupling 
and uncoupling double well potentials. In such double well 
potentials, one may introduce an alternating bias between the 
chemical potentials of neighboring sites. This allows for the 
following three steps: 

(I) Periodic patterns of atoms can now be prepared by iso- 
lating double wells and introducing a bias between odd and 
even sites. Further experimental techniques make sure that 
multiply occupancies are highly suppressed, leaving, e.g., a 
sequence of empty and singly-occupied sites. 

(II) Period-2 local density measurements can be performed 
by mapping odd and even sites to different Brillouin zones: 
Each part of a decoupled double well has multiple bands sep- 
arated by well-defined energies. Biasing the odd sites rela- 
tive to the even ones by an energy in excess of the separation 
energy of the band-separation energy, odd-site particles are 
reloaded into the higher bands of the even sites, whereas the 
even-site particles stay in the lowest band. A time-of-flight 
mapping then reveals the even-site particles in the first, the 
odd-site particles in the higher Brillouin zones. 

(III) Correlations between even and odd sites can be mea- 
sured using more sophisticated techniques lfl6ll . 

Setting and initial condition. - We propose to start from 
a two-periodic initial state prepared by the superlattice setup, 
where the odd sites are occupied by exactly a single boson, 
all even sites being empty, such that the initial state vector 
is |^(0)) = |1, 0,1,0, ...1,0). The 2a-lattice is then sud- 
denly switched off, generating a quenched, non-equilibrium 
situation, and the state vector will evolve in time \ip(t)} = 
e~ ltH / h \ip(0)) according to the Bose-Hubbard Hamiltonian 

L jj L L 

i— 1 i— 1 i— 1 

U and J being the interaction and hopping parameters of the 
Bose-Hubbard model that can be calculated from the lattice 
parameters (J = th = 1). L defines the system size, which 
will be taken to be even. Occupation of higher order bands 
will not be considered, and we stay within the limit of appli- 
cability of the Bose-Hubbard model. We will now see how lo- 
cal relaxation — accompanied by homogenization of densities 
and entanglement — manifests itself in such a periodic setting. 

Signatures of local relaxation for observable quantities. - 
We find that quantitative signatures of local relaxation can be 
obserbed by the global measurement of the total occupation 
of even and odd sites (II), (N 0ye (t)), and correlations between 
nearest neighbors (III), (bi(i)bi+i (t)). In a translationally in- 
variant setting, such measurements amount to accessing lo- 
cal observables as (N e (t)) — 2j=i v*2t(i)) = L(fi2i(t)) /2. 




FIG. 2: Local density (hi(t)) vs. time, showing local relaxation. 
Shown is the time evolution of an even and an odd site for U = and 
U — 1.5, (7 = 3, and U — 8, as well as the estimated negative expo- 
nents of asymptotic power-law decay (right lower figure). Note the 
strong deviation from the non-interacting limit and the strong sup- 
pression of density oscillations for U = 3, as well as the similarity 
for U = 8 with the limit U = oo. 



The limiting cases U = and U = oo are or can be 
mapped to free models. For U = 0, one finds in the Heisen- 
berg picture S<(t) = Y^=i v i-j(t)bj( Q )' where Vj(i) = 

Efe=i e " itAfce27dfeZ/L / LandA fc = -M-2cos(27rfc/L). In the 
limit U — > oo, the interaction manifests itself in that bosons 
become hardcore, rendering them solvable via the Jordan- 
Wigner transformation to spinless fermions. 

To study the relaxation dynamics for finite U, we turn to the 
time-dependent variant |Qj] of DMRG OH, allowing to fol- 
low the coherent time-evolution of strongly interacting quan- 
tum systems very precisely, keeping up to 5000 states in the 
matrix-product state calculations. 

In all of the considered cases, we do find local relaxation: 
The intuition is that the incommensurate influences of trav- 
elling excitations J2, 01 lead to a mixing in time, and hence 
the emergence of properties that locally appear like ones of 
maximum entropy states. These excitations proparate at most 
with the Lieb-Robinson velocity, giving rise to an approxi- 
mate locality in the lattice lfl5ll . The situation is especially 
clear for [7 = 0, where relaxation is due to dephasing in 
the sense that freely propagating excitations lead to reduced 
state contributions of quickly oscillating phases that average 
out J2l ■ One can indeed prove that for any subblock I of con- 
secutive sites, the system relaxes to the reduction of a max- 
imum entropy state given the constraints of motion. This is 
true without time average, arbitrarily exactly to any small er- 
ror in trace-norm, and for an arbitrarily long time lfl9ll . For 
interacting cases with U > 0, we would expect that for very 
short times, observables evolve as in the U = limit, with a 
crossover in behavior for longer time scales when they start 
interacting. We will explore to what interaction strengths the 
remarkably simple limiting pictures remain essentially valid 
and in what regimes one can identify genuinely different re- 
laxation dynamics. 
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FIG. 3: Real and imaginary parts of the correlations to neighbors 
(bi+i(t)bi(t)) as a function of time, for different values of U. The 
inset shows the equilibrated value of the real part of correlations to 
neighbors for large t for different values of U (with a solid line pro- 
portional to 1/U as a guide to the eye). 



Tune evolution of densities. - Signatures of non- 
equilibrium relaxation dynamics are specifically apparent in 
the observation of local densities. Local densities evolve in 
both exactly solvable cases U — 0, oo as 
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for L — > oo, where J n denote the Bessel functions of first 
kind. Odd- and even-site densities relax symmetrically about 
the n = 1/2 axis to n = 1/2, with an asymptotic decay in 
time as t^ 1 / 2 + o(i _1 / 2 ), see Fig. 13] 

For the interacting cases, all t-DMRG results are compati- 
ble with a relaxation of densities to n = 1/2. As excitations 
are responsible for local relaxation, on very short time scales 
(t < 1) particles have typically not interacted yet and are not 
quite sensitive to different values of U, giving rise to a similar 
behavior as for U = 0, oo. The dynamics of relaxation de- 
viates quite strongly for intermediate times, however. This is 
clearly exhibited in Fig. |2] (for L = 32). The decay behavior 
as estimated from the data is also shown in by power-law ex- 
ponents compatible with the data. We clearly see that for very 
small U and U > 4, the limiting non-interacting cases provide 
a very good approximation to the interacting dynamics (one 
encounters a slope similar to —1/2, as for the limiting cases). 
For intermediate values of U, close to the critical point, scat- 
tering appears to be most effective, leading to a strongly en- 
hanced damping and relaxation. Such strong deviations from 
the limiting behavior are expected be visible experimentally, 
as a signature of interaction effects. 

Time evolution of correlators. - We now turn to nearest- 
neighbor correlator (b i+1 (t)bi(t)) . This quantity is specifi- 
cally interesting, going beyond local densities: The build-up 
of correlations in time starting from the uncorrected initial 
state becomes apparent. In the limiting free cases, we find 
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The real part of the correlator is zero for all times, whereas 
the imaginary part relaxes to with an asymptotics of i _1//2 , 
following a quick growth to a maximal value of about 0.28 



at time t ~ 1/2, reflecting the buildup of correlations due to 
particle motion with speed approximately linear in J. 

Fig Sshows the real and imaginary parts of the correlators. 
The buildup of the imaginary part of correlations is largely 
independent of U, reflecting the fact that over the distance 1 
between particles at t = 0, few collisions have yet happened. 
When the interaction becoming visible, the relaxation dynam- 
ics follows quite different paths, local relaxation again being 
fastest around U ~ 3. 

For all finite U, the real part converges to a finite value. 
Indeed, for large U the converged value is well-approximated 
by a U^ 1 curve (for U > 4). This dependence is, in fact, 
exactly what one would expect in the thermal or Gibbs state 
of the Bose-Hubbard model. This can already be seen using 
thermal perturbation theory, from which one obtains 

(*li*> = ^E e "^" +£m) -(-+ 1 )^? if. 

T U z -—' z^yn — vn—l) 

n,m v J 

up to o(J/U), where z = J2n e ~ l3En anc ^ E n — Un(n — 
l)/2 — \m are the local energies of the unperturbed Hamil- 
tonian. So, indeed, within the validity of perturbation theory, 
we do find the anticipated linear dependence on 1/U, as seen 
also in DMRG simulations. This insight further corroborates 
the intuition that locally, the state is indistinguishable from the 
situation as if the system was, globally, in a state maximizing 
the entropy, respecting the constraints of motion J2, Si • 

Entanglement dynamics. - The travelling excitations give 
rise to entanglement between any subsystem and the rest of 
the lattice lflu[l2ll ; indeed to maximal entanglement given the 
constraints of motion if they are local. The speed of infor- 
mation transfer also governs the dynamics at which bipartite 
and long-range entanglement is being built up. The buildup 
of entanglement is both a resource (being ultimately respon- 
sible for local relaxation) and a burden: Linear entanglement 
growth as found here leads to an exponential growth in the 
numerical resources, limiting simulations to t ~ 6 J. 

Time scales of local relaxation. - In all quantities studied 
so far, one can identify three regimes in time: 

(a) Initially, correlations are being built up. In this regime, 
Jt < 1, the dynamics is largely independent of U, as colli- 
sions have not yet become important. 

(b) The second time regime is the one of actual local re- 
laxation. The fast oscillatory dynamics between neighboring 
sites is accompanied by slow local relaxation, resulting from 
the incoming excitations from farther and farther sites, broad- 
ened by dispersion. This results in relaxation, not due to deco- 
herence but due to the dilution of information over the lattice. 
This information propagation is expected to happen at a finite 
speed related to J IU5I1 . in a "ballistic transport". In the free 
models, one finds a polynomial decay: The Bessel function 
fulfills Jo(x) = x" 1 / 2 + o(x~ 1 / 2 ). For finite U, numerics is 
consistent with polynomial decay, which for intermediate U 
seems to be much faster. 

(c) The third and very large time, not yet visible on the 
simulation time scales, is the recurrence time where the finite 
size of the quantum system becomes visible. 

The interaction strength also marks three regimes: 
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(A) For small interactions U < 1 the dynamics strongly 
resembles the non-interacting bosonic limit. 

(B) For large values of U the dynamics is very similar to 
the hardcore bosonic or free fermionic limit, with relaxation 
exponents being similar to —1/2. The observed local corre- 
lations proportional to 1/U are consistent with assuming that 
we locally observe a global state having maximum entropy. 

(C) The case of intermediate U ~ 3 appears to mark the 
crossover between the two free cases of U — and U = oo, 
characterized by the most efficient relaxation. 

Quasi-momentum distribution. - Remarkably, the quasi- 
momentum distribution (QMD), measurable via time-of- 
flight, defined as S(q, t) = E^=i e^"^ {b\(t)bj (*)) /L can 
be shown not to relax for U = 0. This creates, for small U, 
the interesting situation that local relaxation can be probed, 
while at the same time signatures of the memory of the initial 
condition can be certified by the absence of relaxation for the 
QMD j20h . 

Effect of a harmonic confinement potential. - We find little 
influence on the local quantities, up to the finite-size recur- 
rence time, which can be shortened in the presence of a re- 



alistic trap: the generated excitations no longer travel with a 
constant speed, but are slowly reflected. 

Summary. - In this work, we have introduced a setting 
allowing for the observation of apparent thermalization in 
closed quantum many-body systems, without the need of ad- 
dressing single sites of a lattice system. We show how, making 
use of optical superlattices, signatures of local relaxation, in- 
cluding relevant time scales and dependencies on interactions, 
can be probed. In such dynamics, local maximum entropy 
states — Gibbs states in the absence of further constraints — 
emerge without having any external heat bath: Instead, in a 
symmetric fashion, each system forms the effective environ- 
ment of the other in complex many-body dynamics. It is the 
hope that ideas along the line of the present work open up a 
way of experimentally quantitatively simulating an instance 
of the fundamental physical process of thermalization. 
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